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Over the last years a new generation of model atmosphere codes, which include the effects of 
, metal line-blanketing of millions of spectral lines in NLTE, has been used to re-determine the 

' properties of massive stars through quantitative spectral analysis methods applied to optical, 

, IR and UV spectra. This has resulted in a significant change of the effective temperature scale 

04 ■ of early type stars and a revision of mass-loss rates. Observed mass-loss rates and efi'ective 

temperatures depend strongly on metallicity, both in agreement with theoretical predictions. The 
new model atmospheres in conjunction with the new generation of lOm-class telescopes equipped 
with efficient multi-object spectrographs have made it possible to study blue supergiants in 
' galaxies far beyond the Local Group in spectroscopic detail to determine accurate chemical 

composition, extinction and distances. A new distance determination method, the flux weighted 
gravity - luminosity relationship, is discussed as a very promising complement to existing stellar 
distance indicators. 

Observationally, there are still fundamental uncertainties in the determination of stellar mass- 
loss rates, which are caused by the fact that there is evidence that the winds are inhomogeneous 
, and clumped. This may lead to major revisions of the observed rates of mass-loss. 

o 

O ■ 1. Introduction 

■ Hot massive stars are cosmic engines of fundamental importance not only in the local 
Oh| but also in the early universe. A first generation of very massive stars has very likely 
Q ■ influenced the formation and evolution of the first building blocks of galaxies. The spectral 
| H I appearance of Lyman break galaxies and Lya-emitters at high redshift is dominated by an 
^ ■ intrinsic population of hot massive stars. Gamma-ray bursters are very likely the result 

of terminal collapses of very massive stars and may allow to trace the star formation 
history of the universe to extreme redshifts. 

It is obvious that the understanding of important processes of star and galaxy for- 

■ mation in the early universe is intimately linked to our understanding of the physics of 
I massive stars. The observational constraints of the latter are provided by quantitative 

spectroscopic diagnostics of the population of hot massive stars in the local universe. It is 
the goal of this contribution to provide an overview about the dramatic progress, which 
has been made in this field over the last five years. 

There are two factors which have contributed to this progress, new observational fa- 
cilities such as the optical/infrared telescopes of the lOm-class on the ground and obser- 
vatories in space allowing for spectroscopy in the UV (HST, FUSE, GALLEX), IR (ISO, 
Spitzer) and at X-ray wavelengths (XMM, Chandra) and the enormous advancement of 
model atmosphere and radiative transfer techniques. As for the latter, it is important 
to realize that modelling the atmospheres of hot stars is a tremendous challenge. Their 
physics are complex and very different from standard stellar atmosphere models. They 
are dominated by a radiation field with energy densities larger than or of the same order 
as the energy density of the atmospheric matter. This has two important consequences. 
First, severe departures from Local Thermodynamic Equilibrium of the level popula- 
tions in the entire atmosphere are induced, because radiative transitions between ionic 
energy levels become much more important than inelastic collisions. Second, supersonic 
hydrodynamic outflow of atmospheric matter is initiated by line absorption of photons 
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Figure 1 . Evolutionary tracks of massive stars in the HRD with (sohd) and without (dashed) 
stellar rotation, [Meynet fc Maeder (2000)| ) and the domains of spectral types discussed in this 
review. 



transferring outwardly directed momentum to the atmospheric plasma. This latter effect 
is responsible for the existence of the strong stellar winds observed and requires the use 
of NLTE model atmospheres, which include the hydrodynamic effects of stellar winds 
(see Kudritzki (1998)_ for a detailed description of the physics of hot star atmospheres). 

The winds of hot massive stars are fundamentally important. Their energy and momen- 
tum input into the ISM is significant creating circumstellar shells, stellar wind bubbles 
and initiating further star formation. They affect the stellar evolution by modifying evo- 
lutionary timescales, chemical profiles, surface abundances and stellar luminosities. They 
also have substantial effects on the structure of the stellar atmospheres. They dominate 
the density stratification and the radiative transfer through their transonic velocity fields 
and they modify the amount of the emergent ionizing radiation significantly (Gabler et 
al., 1989, 1991, 1992, Najarro et al., 1996). 

While the basic concepts for the hydrodynamic atmospheres of hot stars and the 
spectroscopic diagnostics of their parameters and stellar winds have been developed in 



the eighties and nineties (see reviews by Kudritzki & Hummer (1990) Kudritzki (1998) 



and Kudritzki & Puis (2000) I, the development of the most recent generation of model 
atmospheres, which for the first time accounts selfconsistently for the effects of NLTE 
metal line-blanketing (see section 2), has lead to a dramatic change of the diagnostic 
results. A new effective temperature scale has been obtained for the 0-star spectral types, 
which are apparently significantly cooler than originally assumed and, thus, on average 
have a lower luminosity and mass and also provide less ionizing photons, than previously 
thought. We will discuss these effects in detail in section 2, 3, and 4. Note that we will 
focus our review on "normal" hot massive stars in a mass-range between 15 and 100 Mq 
in well-established evolutionary stages such as dwarfs, giants, and supergiants (Fig.^. 
We will not discuss objects with extreme winds such as Wolf-Rayet stars, luminous blue 
variables in outbursts, etc. For those, we refer the reader to the contributions by Paul 
Crowther and Nathan Smith in this volume. 

The new model atmospheres in conjunction with the new generation of lOm-class tele- 
scopes equipped with efficient multi-object spectrographs allow to study blue supergiants 
in galaxies far beyond the Local Group in great spectroscopic detail to determine accu- 
rate chemical composition, extinction and distances. We will report recent results in this 
new field of "extragalactic stellar astronomy" in section 5. In section 6, we will discuss a 
new distance determination method based on stellar photospheric spectroscopy, the flux 
weighted gravity - luminosity relationship. 
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Figure 2. Model UV spectra of an integrated population of massive stars in a starburst 
gal axy displayed as a function of metallicit y. The model atmosphere code used was WMBa- 
sic ( [Pauldrach , Lennon fc Hoffmann (2001)*) and the population synthesis was done with Star- 
burst99. (1^'rom Rix, Pcttini, Lcithcrcr ct al. (2004) J . 

The diagnostics of stellar winds, in particular of mass-loss rates, are also affected by the 
proper accounting for line-blanketing affects. While the general scaling relations of stellar 
wind parameters with stellar luminosity, mass, radius, and metallicity as described in 
Kudritzki fc Puis (2000) remain qualitatively unchanged, important quantitative changes 



have been found. There are still puzzling uncertainties with regard to the observed rates 
of mass-loss related to the inhomogeneous structure of the stellar wind outflows. We will 
present and discuss the most recent results in section 7. 

As indicated above, the new generation of model atmosphere codes has important ap- 
plications on the interpretation of spectra of star-bursting galaxies in the early universe. 



Fig.[21is a nice example. This work by Rix, Pettini, Leitherer et al. (2004) has been used 
to constrain star formation rates and metallicities in high redshift Lyman-break galax- 
ies. Another example is the work by Barton et al. (2004) which uses stellar atmosphere 
model predictions by Bromm, Kudritzki & Loeb (2001) and 'Schaerer (2003) to explore 
the possibility to detect the first generation of very massive stars at redshifts around 
ten with present day 8m-class telescopes and with future (diffraction limited) optical/IR 
telescopes of 30m aperture such as the GSMT (see also report by the GSMT Science 
Working Group, Kudritzki et al. (2003a) I. As it turns out, not only would the GSMT be 
able to detect such objects it would also allow to constrain the IMF of massive stars in 
the early universe from the relative comparison of and HeII1640 recombination lines 
caused by the population of the first stars. 



2. The effects of metal line blanketing 

The inclusion of the opacity of millions of spectral lines in NLTE has two major effects. 
First, it changes the spectral energy distribution in the UV because of strong metal line 
absorption in the outer atmosphere ("line-blanketing"). An example is given in Fig. |2| 
While the UV flux is reduced, the physical requirement of radiative equilibrium and 
conservation of the total flux leads to an increase of the emergent flux at optical and 
IR wavelengths relative to the unblanketed case. This increase comes from the fact that 
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Figure 3. Effects of metal line blanketing. Left: Emergent flux of a NLTE hydrodynamic 
0-star model of Te//=40,000K with (solid) and without (dashed) metal line blanketing. Right: 
Local kinetic temperature as a function of monochromatic optical depth at blue wavelength for 
the same two models. The photospheric backwarming as discussed in t he text as the result of 
blanketing is clearly seen. The NLTE code FASTWIND ( |Puls et al. "(2005)) was used for the 
calculations. 



a significant fraction of the photons absorbed in the outer layers are emitted back to 
the inner photosphere providing additional energy input and, thus, heating the deeper 
photospheres. This second "backwarming" effect increases the local kinetic temperature 
(see Fig.|2)l, which then leads to stronger photospheric continuum emission and modifies 
diagnostically important ionization equilibria such as Hel/II, which are used for the de- 
termination of Teff- Fig. ^demonstrates how the Hel/II ionization equilibrium is shifted 
towards lower Teff because of the backwarming effect. At the same time, the pressure- 
broadened wings of the Balmer lines (the standard diagnostic for log g) become weaker, 
because the millions of metal lines increase the radiative acceleration grad and decrease 
the effective gravity geff = g - grad- As a result higher gravities are needed to fit the 
Balmer lines (see Fig. 0J in addition to the lower temperatures obtained from the he- 
lium ionization equilibrium. In summary, the use of blanketed models leads to systematic 
shifts in the (log g, log Teff) - plane towards lower temperatures and higher gravities. For 
stars with known distances to be analyzed this means lower luminosities and lower masses. 
We also note that the presence of dense stellar wind envelopes increases the effects of 
backwarming (see Hummer (1982)| Abbott fc Hummer (1985)) . In addition, stellar winds 
also affect the strengths of diagnostically crucial Hel absorption lines by contaminating 
the absorption with stellar wind line emission ( Gabler, Gabler, Kudritzki et al. (1989)[ 



Sellmaier et al. (1993) Repolust et al. (2004) K 



The combined effects of line blanketing and backwarming affect also the ionizing 
fluxes. Amazingly, for the ionization of hydrogen the changes are very small as the ef- 
fects of blanketing and backwarming balance each other. However, the ionization of ions 
with absorption edges shorter than the one for hydrogen is significantly affected (see 
Kudritzki (2002)[ [Martins et al. (2005a) 1. 



Over the last years very powerful and user friendly software packages to calculate NLTE 
model atmospheres with metal line blanketing have been developed, which are now avail- 
able to the astronomical community and which have been intensively applied for spectral 
diagnostics of hot massive stars. At Munich University Observatory the codes "WM- 
Basic" ( |Pauldrach, Lennon fc Hoffmann (2001)] ) and "FASTWIND" ( |Puls et al. (2005)| | 
were developed in partial collaboration with the authors of this article and at Pitts- 
burgh the code "CMFGEN" is the result of intensive work by John Hillier and coUabo- 
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Figure 4. Effects of metal line blanketing on the stellar diagnostics. Plotted are the isocontours 
of log [WA(HeI4471)/WA(HeII4542)] (left) and Wa(Ht.) (right) in the (log g, log Te//) - plane. 
Dashed isocontours are unblanketed models, solid are blanketed. Vectors indicate the shifts 
caused by the effects of NLTE metal line blan keting. As for the p revious figures the calculations 
were done with the NLTE code FASTWIND | |Puls et al. (2005)| l. 
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Figure 5. Effects of metal line blanketing on the stellar diagnostics. Left: Shifts of the location 
of individual O-stars in the (log g, log Te//)-plane caused by the use of metal line blanketed 
models in the analysis of the hydrogen and helium optical line spectrum. Righ t: Shifts of modified 
stellar wind momenta (see section 7) as a function of luminosity. (From ^Repolust et al. (2004)| l. 



rators ( Hillier et aL (2003) I . Ivan Hubeny and Thierry Lanz have developed "TLUSTY" 
iLanz & Hubeny (2003) I, which does not include the hydrodynamics of winds but can 
be apphed in all cases, where the stellar winds are weak. 



3. The effective temperature scale of O-stars 

With the new generation of NLTE metal line blanketed model atmospheres available 
quite a number of very detailed spectroscopic studies of O-stars have been published 
over the last five years. While the diagnostic methods used are generally still the same as 
in the earlier work by [Kudritzki (1980)| [Kudritzki et al. (1983)[ [Kudritzki et al. (1989)| 



Kudritzki &: Hummer (1990)[ Herrero, Kudritzki, et al. (1992)[ which was based on a 



first generation of hydrostatic NLTE model atmospheres with hydrogen and helium opac- 
ity only, the results coming out of the new work are substantially different because of the 
tremendous improvements of the model atmospheres. Since the effects of metal line blan- 
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Figure 6. Left: Effective temperature as a function of spectral type for Milky Way 0-stars. 
Supergiants are plotted as circles and giants and dwarfs as squares. The new scales are rep- 
resented by the solid curves, and the old scales by the dashed curves. For discussion see text. 
Right: The new effect ive temperat ure scale of O-stars in the SMC (solid) relative to the Milky 
Way (dashed). (From |Massey, Puis, Pauldrach et al. (2005)} . 



keting depend on the stellar metallicity, we will discuss 0-stars in the solar neighborhood 
of the Milky Way and in the Magellanic Clouds separately. 



3.1. 0-stars in the Milky Way 
Milky Way 0-stars have been studied by|Pauldrach, Lennon & Hoffmann (2001) Herrero et al. (2002)| 



Martins et al. (200 2) , Bia nchi fc Ga rcia (20 02)[|Garcia fc Bian chi (200 4)||Repolust et al. (2004)[ 



Markova et al. (2004)nMartins et al. (2005a)HMokiem et al. (2005)nBouret et al. (2005) 



The result of all these studies with regard to the effective temperature scale is quite dra- 
matic. An example is given by Fig. El which displays the shifts in the (log g, log Tg//)- 
plane of indivual 0-stars, when the classical analysis of hydrogen and helium lines based 
on unblanketed model atmospheres is replaced by blanketed models. As to be expected 
from the previous section, the use of blanketed models leads to cooler effective temper- 
atures, which is simply caused by the fact that blanketed models have higher intrinsic 
local photospheric temperatures and thus models with lower Te/ / are required to fit the 
observed helium ionization equilibrium. In addition, blanketed models have a lower local 
gas pressure and, thus, a higher gravity is generally needed to fit the pressure broadened 
Balmer lines. 

While systematic effects towards somewhat lower effective temperatures and higher 
gravities, were always expected in previous work based on unblanketed models, the large 
shifts as displayed in Fig. Income as a surprise. We note that the effects are strongest for 
very hot objects and low gravity supergiants, which have very strong and dense winds af- 
fecting the helium ionization equilibrium additionally by the mechanisms discussed above. 
In some cases we encounter effective temperature changes of the order of 10 to 15 per- 
cent. Based on the work by Repolust et al. (2004)[ Massey, Puis, Pauldrach et al. (2005) 



have introduced a new effective temperature scale for Milky Way 0-stars, which is dis- 
played in Fig. Eland compared to the old scale by Vacca et al. (1996) We realize that, 
in particular for supergiants, the differences between the two scales are dramatic. A 
very similar new effective temperature scale has also been introduced independently by 



Martins et al. (2005a) based on a comprehensive study of a large sample of 0-stars. 
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Figure 7. Model atmosphere fits of IR hydrogen and hehum lines to determine stellar param- 
eters. Note that Hel 2.11/i is blended by nitrogen lin es not included in the c alculations. Br^, is 
blended by HII region emission in some cases. (From [Repolust et al. (2005)| l . 

3.2. Quantitative IR spectroscopy of 0-stars 

Massive stars are frequently found in star forming regions heavily obscured by interstel- 
lar dust. In such a situation IR spectroscopy is the only way to obtain information. 
In the galactic cente r [Najarro et al. (1994)[ [Najarro et al. (1997)[ |Figer et al. (1998)| 
Najarro et al. (2004) (see also the contribution by Don Figer in this volume) have demon- 
strated how quantitative IR spectroscopy can be used to determine the stellar parameters 
of very extreme supergiants and Wolf-Rayet stars. Similar work has now been carried 
more recently for "normal" 0-stars with the goal to find out whether the analysis of 
the hydrogen and helium lines in the IR yields accurate information about stellar pa- 
rameters consistent with the ones obtained from the analysis of the optical spectrum. 
Hanson, Kudritzki, Kenworthy et al. (2005) and Repolust et al. (2005) have used Sub- 



aru H and K band spectra of high S/N to study a large sample of 0-stars with very 
encouraging results showing that IR spectroscopy leads to effective temperatures, gravi- 
ties and mass-loss rates practically identical with those determined from optical spectra 



(Fig. [7| gives an example). In parallel, Lenorzer et al. (2004) have carried out a system- 
atic model atmosphere study of the IR spectral diagnostics for normal 0-stars focussing 
on both stellar lines as well as HII region emission lines. It is obvious from this new work 
that for future investigations of massive stars in dense and highly obscured star forming 
regions quantitative IR spectroscopy has an enormous potential. 

3.3. Effects of metallicity. The effective temperature scale of 0-stars in the Magellanic 

Clouds 



From the discussion in section 2 it is clear that metallicity must have an influence on the 
strengths of blanketing effects and, thus, the effective temperature scale of 0-stars. The 
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Figure 8. Model atmosphere fits of optical h ydroRen and helium lines of an O-star in t he 
LMC to determine stellar parameters. (From |Massey, Bresolin, Kudritzki et al. (2004)| l. 

ideal laboratory to investigate the metallicity dependence are the Magellanic Clouds be- 
cause of their lower metallicity relative to the Milky Way. jMassey, Bresolin, Kudritzki et al. (2004) 



and Massey, Puis, Pauldrach et al. (2005) and, independently, Mokiem et al. (2004)| and 



Mokiem et al. (2006a) have carried out a systematic and comprehensive spectroscopic 
study of 0-stars in the clouds. Fig. El and Fig. |S1 describe the work done by Massey and 
collaborators. From a detailed fit of the optical hydrogen and helium lines and HST 
and FUSE UV spectra they were able to determine stellar parameters and stellar wind 
properties. The result with regard to the effective temperature scale is striking and in 
agreement with the expectation based on model atmosphere theory. For a given spectral 
type defined by the relative strengths of Hel and Hell lines 0-stars at lower metallicity 
are significantly hotter than their galactic counterparts, a result which is also clearly con- 
firmed by the work of Mokiem and collaborators. For the determination of IMFs in other 
galaxies based on spectral classification or for the investigation of the stellar content of 
galaxies or Hll-regions based on nebular emission line analysis this is an important effect 
to be taken into account. Also for the study of integrated spectra of starburst galaxies 
the metallicity dependence of the effective temperature scale might be important. 

3.4. Systematic effects depending on the analysis method 

In the work discussed so far effective temperatures and gravities have been determined 
from the fit of the He I/II ionization equilibrium of optical or IR helium lines and the 
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Figure 9^ O-star effective temperature scales according to 

[Massey, Puis, Pauldrach et al. (2005) | for the Milky Way (left) and the the SMC (right) 
based on the analysis of the Hel/11 ionization equilibrium at optical wavelengths. Overplotted 
are the results of spectral analyses using UV metal lines as temperature indicator (see text for 
the data sources). Supergiants are plotted as circles giants and dwarfs as squares. 



Balmer lines. Since 0-stars are UV bright and hundreds of photospheric lines can be 
identified in high resolution UV spectra (lUE, HST, ORFEUS, FUSE), it is tempting 
to use this spectral range to obtain independent information about the stellar parame- 
ters. As it turns out, it is difficult to constrain gravities solely through the UV, however, 
there are ionization equlibria such as Fe IV/V/VI or CIII/IV which can be used for a 
determination of temperatures. For Milky Way 0-stars such work has been carried out by 
Pauldrach, Lennon fc Hoffmann (2001)| IBianchi fc Garcia (2002)|[Garcia fc Bianchi (2004) 



and Bouret et al. (2005) As shown in Fig.|51significantly lower Tg// are obtained in some 
cases from the UV lines compared to optical/IR Hel/II lines. In particular, the work by 
Bianchi & Garcia (2002) and Garcia & Bianchi (2004) has resulted in very low effective 



temperatures. 

A similar result has been obtained by Heap, Lanz & Hubeny (2006) for SMC 0-stars 
mostly of luminosity class V. Fig.|5|shows their results compared to the effective tempera- 



ture scale obtained by Massey, Bresolin, Kudritzki et al. (2004) and Massey, Puis, Pauldrach et al. (2005) 
Fig. El also displays the results found by [Crowther et al. (2002)[ [Hillier et al. (2003)[ 



Bouret et al. (2003)[ and Evans et al. (2004)[ which show a less extreme but similar 



trend. 

This indication of a systematic effect depending on the analysis method deserves a 
very careful and systematic future investigation. What is needed is a comprehensive 
simultaneous UV, optical, IR analysis of all the 0-stars studied so far. The observational 
material seems to be available or can be easily obtained. Of course, this will be a time 
consuming effort requiring a lot of detailed spectroscopic and model atmosphere work. 
However, it is extremely important that the reasons for these obvious discrepancies are 
well understood. 



4. The effective temperature scale of B and A supergiants 

As shown in Fig. ^ massive 0-stars evolve into B and later A supergiants. Based on 
unblanketed, hydrostatic NLTE model atmospheres an effective temperature scale was 
introduced for galactic B supergiants by McErlean et al. (1998, 1999). More recently. 
Trundle et al. (2004)[|TVundle fc Lennon (2005)1 and lCrowther, Lennon fc Walborn(2005) 



have analyzed a large sample of Milky Way and SMC B supergiants with the improved 
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Figure 10. The effective temperature scales of early B-supergiants (left, the dashed curve shows 
the old scale based on unblanketed NLTE models) and late B and early A-supergiants (right) 
of luminosity class la. Note that in the right plot on the abcissa the value of -2 corresponds to 
spectral type B8. 



models and have revised the effective temperature scale for the earliest spectral types of 
luminosity class la. The results are shown in Fig.^| Note that if the spectral classification 
scheme introduced by |Lennon, D.J. (1997) is used, the effective temperature scale should 
not depend on metallicity. This is indeed the case, as a comparison of the SMC sample 
I Trundle k Lennon (2005]| ) with the MW sample ( |Crowther, Lennon fc Walborn(2065) I 
shows. 



Kudritzki, Bresolin, & Przybilla (2003) have recently introduced an effective tempera- 



ture scale for late B and early A supergiants of Iximinosity class la and of solar metallicity, 
which is also displayed in Fig.^] A discussion of the metallicity dependence can be found 



Evans fc Howarth (2003) 



5. Extragalactic stellar astronomy with B and A supergiants 

As shown in Fig. ^ 0-stars with masses between 15 to 40 Mq evolve into super- 
giants of spectral B and A. Although not as massive and luminous as the most massive 
0-stars, these stars are the brightest "normal" stars at visual light with absolute mag- 



nitudes —7.0 ^ My ^ —9.5, see Bresolin (2003) (By "normal" we mean stars evolving 



peacefully without showing signs of eruptions or explosions, which are difficult to han- 
dle theoretically and observationally) . While 0-stars emit most of their radiation in the 
extreme and far UV because of their high atmospheric temperature, late B and early 
A supergiants are cooler and their bolometric corrections are much smaller because of 
Wien's law so that their brightness at visual light reaches a maximum value during stel- 
lar evolution. It is this enormous intrinsic brightness at visual light, which makes them 
extremely interesting for extragalactic studies far beyond the Local Group. 

During their smooth evolution from the left to the right in the HRD massive stars 
are crossing the temperature range of late B and early A-supergiants in a timescale 
on the order of several 10'^ years (Meynet & Maeder (2000) I. During this short evo- 
lutionary phase stellar winds with mass- loss rates of the order 10~^ M0yr~^ or less 
( Kudritzki & Puis (2000) I do not have enough time to reduce the mass of the star sig- 
nificantly so that the mass remains constant. In addition, as Fig.^shows, the luminosity 
stays constant as well. The fact that the evolution of these objects can very simply be 
described by constant mass, luminosity and a straightforward mass-luminosity relation- 
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ship makes them a very attractive stehar distance indicator, as we will explain later in 
this review. 

As evolved objects the blue supergiants are older than their 0-star progenitors, with 
ages between 0.5 to 1.3 x 10^ years (Meynet & Maeder (2000)). All galaxies with ongoing 
star formation or bursts of this age will show such a population. Because of their age 
they are spatially less concentrated around their place of birth than 0-stars and can 
frequently be found as isolated field stars. This together with their intrinsic brightness 
makes them less vulnerable as distance indicators against the effects of crowding even 
at larger distances, where less luminous objects such as Cepheids and RRLyrae start to 
have problems. 

With regard to the crowding problem we also note that the short evolutionary time of 
10^ years makes it generally very unlikely that an unresolved blend of two supergiants 
with very similar spectral types is observed. On the other hand, since we are dealing 
with spectroscopic distance indicators, any contribution of unresolved additional objects 
of different spectral type is detected immediately, as soon as it affects the total magnitude 
significantly. 

Thus, it is very obvious that blue supergiants seem to be ideal to investigate the 
properties of young populations in galaxies. They can be used to study reddening laws 
and extinction, detailed chemical composition, i.e. not only abundance patterns but also 
gradients of abundance patterns as a function of galactocentric distance, the properties 
of stellar winds as function of chemical composition and the evolution of stars in different 
galactic environment. Most importantly, as we will demonstrate below, they are excellent 
distance indicators. 



6. Quantitative stellar spectroscopy beyond the Local Group 

Enormous progress has been made over the last years in the development of accurate 
NLTE spectral diagnostics of B and A supergiants. Using high resolution and high S/N 
spe ctra|Urbaneja (2004), Urbaneja et al. (2005b)||Trundle etal. (2004),,Trundle fc Lennon (2005)1 



and Crowther, Lennon & Walborn(2005) have studied early B-supergiants in Local Group 
galaxies and the Milky Way to determine stellar properties, chemical composition and 
abundance gradients. For late B and early A supergiants Przybilla et al. (2001a, b,c, 
2006) have developed very detailed model atoms for the NLTE radiative transfer diag- 
nostics, which allow for an extremely accurate determination of effective temperature 
(1%), gravities (0.05 dex), and chemical abundances (0.1 dex). Detailed abundance stud- 
ies of A supergiants in many Local Group galaxies using these NLTE methods were 
carried out by Venn et al. (1999, 2000, 2001, 2003) and |Kaufer et al. (2004")] The anal- 
ysis technique is similar to 0-stars, except that different ionization equilibria are used 
for the determination of effective temperatures (Sill/III/IV for early B supergiants and 



OI/II, NI/II, Mgl/II, SII/III for late B and early A supergiants, see [Kudritzki (2003)| for 
a more detailed description). 

6.1. Chemical composition 

For extragalactic applications beyond the Local Group spectral resolution becomes an 
issue. The important points are the following. Unlike the case of late type stars, crowding 
and blending of lines is not a severe problem for hot massive stars, as long as we restrict 
our investigation to the visual part of the spectrum. In addition, it is important to 
realize that massive stars have angular momentum, which leads to usually high rotational 
velocities. Even for A-supergiants, which have already expanded their radius considerably 
during their evolution and, thus, have slowed down their rotation, the observed projected 
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Figure 11. Quantititative spectral analysis of an A-supergiant in the Sculptor spiral 
galaxy NGC 300 at a distance of 2 Mpc. The inset shows the fit of the Balme r jump, 
which is used to determine Te//. (From [Kudritzki, Urbaneja, Bresolin et al. (2006b)| see also 
[Bresolin et al. (2002a)l . 



rotational velocities are still on the order of 30kms~^ or higher. This means that the 
intrinsic full half-widths of metal lines are on the order of 1 A. In consequence, for the 
detailed studies of supergiants in the Local Group a resolution of 25,000 sampling a line 
with five data points is ideal. This is indeed the resolution, which has been applied in 
most of the work referred to above. 

However, as we have found out empirically ( Przybilla (2002) I, degrading the resolution 
to 5,000 (FWHM = 1 A) has only a small effect on the accuracy of the diagnostics, as 
long as the S/N remains high (i.e. 50 or better). Even for a resolution of 2,500 (FWHM 
= 2 A) it is still possible to determine Te// to an accuracy of 2 percent, log g to 0.05 dex 
and individual element abundances to 0.1. 

Bresolin et al. (2001, 2002a, 2002b, 2 003 , 2006), Urban eja et a l. (2003, 2005a, 2006), 
Kudritzki, Bresolin, & Przybilla (2003) and Kudritzki, Urbaneja, Bresolin et al. (2006b) 
have used FORS at the VLT with a resolution of 1,000 (FWHM = 5 A) to study blue 
supergiants far beyond the Local Group. The accuracy in the determination of stellar 
properties at this rather low resolution is still remarkable. The effective temperature (for 
late B and A supergiants now determined from the Balmer jump rather than from ioniza- 
tion equilibria) is accurate to roughly 4 percent and the determination of gravity based on 
fitting the broad Balmer lines remains unaffected by the lower resolution and is still good 
to 0.05 dex. Abundances can be determined with an accuracy of 0.2 dex. In Fig. 1111 IT^ 
and 1 131 we show examples of the detailed spectral fits that can be accomplished. Fig. 1131 
demonstrates, how important information about the metallicity gradients of the young 
stellar population in spiral galaxies can be obtained directly from the spectral analysis 
of blue supergiants. 
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Figure 12. Quantititative spectral analysis of an early B-supergiant in the Sculptor spiral 
galaxy NGC 300. (From fUrbaneja et al. (2005a)t . 
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Figure 13. Left: the metallicity gradient of the young stellar population in the Sculp- 
tor spiral gal axy NGC 300 as detemined for quantitative stellar spectroscopy. (From 
[Urbaneja et al . (2005a), Bresolin et al. (2002a, b). Right: Quantitative spectral analysis of an 
A-supergiant in the spiral galaxy NGC 3621 at a distance of 7 Mpc. (From |Bresolin et al. (2001)| l. 



6.2. Extragalactic distance determinations with the "Flux Weighted Gravity - 
Luminosity Relationship (FGLR)" 

The best established stellar distance indicators, Cepheids and RRLyrae, suffer from two 
major problems, extinction and metallicity dependence, both of which are difficult to 
determine for these objects with sufficient precision. Thus, in order to improve distance 
determinations in the local universe and to assess the influence of systematic errors there 
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Fig ure 14. Left: The FG LR of stellar evolution models from [Meynet eFa l. (1994^ 
and [Meynet fc Maeder (2000)] Circles correspond to models with rotation, squares represent 
models without the eitects of rotation. Solid symbols refer to galactic metaUicity and open sym- 
bols represent SMC metaUicity. The solid curve corresponds to a = —3.75 in the FGLR. Note 
that gravities are in cgs units and the temperatures are in units of lO^K. Right: Fit of the higher 
Balmer hues of an A-supergiant in the Sculptor galaxy NGC 300 using two atmospheric models 
with Te// = 9500 K and logg — 1.60 (thick hne) and 1.6 5 (thin line), respectively. The data were 
taken with FORS at the VLT. For further discussion see [Kudritzki, Bresolin, fc Przybilla (2003)| 



is definitely a need for alternative distance indicators, which are at least as accurate 
but are not affected by uncertainties arising from extinction or metaUicity. Blue super- 
giants are ideal objects for this purpose because of their enormous intrinsic brightness, 
which makes them available for accurate quantitative spectroscopic studies even far be- 
yond the Local Group using the new generation of 8m-class telescopes and the extremely 
efficient multi-object spectrographs attached to them (see previous subsection). Quanti- 
tative spectroscopy allows us to determine the stellar parameters and thus the intrinsic 
energy distribution, which can then be used to measure reddening and the extinction 
law. In addition, metaUicity can be derived from the spectra. We emphasize that a reli- 
able spectroscopic distance indicator will always be superior, since an enormous amount 
of additional information comes for free, as soon as one is able to obtain a reasonable 
spectrum. 

A very promising spectroscopic distance determination method based on simple stellar 
physics is the Flux-Weighted Gravity - Luminosity Relationship (FGLR), which was 



introduced by Kudritzki, Bresolin, & Przybilla (2003) When discussing Fig.^in Sect. 1 



we noted that massive stars evolve through the domain of blue supergiants with constant 
luminosity and constant mass. This has a very simple, but very important consequence 
for the relationship of gravity and effective temperature along each evolutionary track. 
From 

L oc R^T^ff = const.; M = const. (6.1) 

follows immediately that 

M ocg ocL (g/T^s) = const. (6.2) 

This means that each object of a certain initial mass on the ZAMS has its specific 
value of the "flux-weighted gravity" g/T^^j during the blue supergiant stage. This value 
is determined by the relationship between stellar mass and luminosity, which to a good 
approximation is a power law 

L (X . (6.3) 



R. P. Kudritzki & M.A. Urbaneja: Parameters and Winds of Hot Massive Stars 15 





-10 
-9 


- • 

+ 

- □ 
o 
* 

o 

A 

r 


NGC300 

NGC68S2 

M31 

Milky Way 

NGC36ai 

LMC 

SMC 
M33 










* 




• 

o 




(mag) 


-8 










• 


• 

• 












-7 
-6 




^ O 


• 

• 


• 


• 


• 

• 










Id 

E 


0.5 













« o 










0.0 
-0.5 




o 








oox>x>o 


» o 


<^ oh 






m 
m 






o 

o 


o 

o 


8- 

o 






o 






o 





Figure 15. The FGLR of B and A supergiants in Local Group galaxies and in the spiral 
galaxies NGC 300 and NGC 3621 at a distance of 2 and 7 Mpc, respectively. The abcissa 
is the same as in Fig. | 14 | (left). (From [Kudritzki, Urbaneja, Bresolin et al. (2006b)| see also 
[Kudritzki, Bresolin, fc Przybilla (2003)| l. J^'or a discussion, see text. 



Inspection of evolutionary calculations with mass-loss, cf. Meynet et al. (1994) and 
Meynet &: Maeder (2000)'{ shows that a; = 3 is a good value in the range of luminosities 



considered, although x changes towards higher masses. With the mass - luminosity power 
law we then obtain 

Li-- (X {glTt^r , (6.4) 
or with the definition of bolomctric magnitude Mboi fx — 2.51ogL 

-Mboi-alog(5/r,tf)+6. (6.5) 

This is the FGLR of blue supergiants. Note that the proportionality constant a is given 
by the exponent of the mass - luminosity power law through 

a = 2.hxl{\-x) . (6.6) 

and a = —3.75 for x = Z. Mass-loss will depend on mctallicity and therefore affect 
the mass - luminosity relation. In addition, stellar rotation through enhanced turbulent 
mixing might be important for this relation. In order to investigate these effects we have 



used the models of Meynet et al. (1994) and Meynet & Maeder (2000) to construct the 
stellar evolution FGLR, which is displayed in J^'ig. 114:1 'i'he result is verv encouraging. 
All different models with or without rotation and with significantly different metallicity 
form a well defined very narrow FGLR. 



In order to verify the existence of the theoretically predicted FGLR Kudritzki, Bresolin, & Przybilla (2003) 
re-analyzed a large sample of late B and early A supergiants in several Local Group galax- 
ies (high resolution spectra) and in the spiral galaxies NGC 300 (2Mpc) and NGC 3621 
(7Mpc) using VLT FORS low resolution spectra as descibed above. Fig. ITU demonstrates 
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how precisely the gravities can be determined even at low resolution. The effective tem- 
peratures were obtained from the spectral types using the relation displayed in Fig. 1101 
Taking into account the valid criticism by Evans & Howarth (2003) about the metallicity 
dependence of the temperature vs. spectral type relationship, Kudritzki, Urbaneja, Bresolin et al. (2006b) 



used the information about the Balmer jump in the VLT/FORS spectrophotometric data 
(see Fig. Illf) to determine effective temperature and metallicity independently. This new 
procedure lead to the FGLR shown in Fig. ^| A least square fit yields a=-3.31 and 
b=13.09 with a one a scatter of 0.24 mag. Fixing the slope to the theoretical value 
a=-3.75 we obtain as a zero point b=-13.73 with (t=0.25. 

The FGLR as displayed in Fig. E| is an extremely tight relationship with a scatter 
comparable to the observed period-luminosity relationships of Cepheids. We conclude 
that blue supergiants provide a great potential as excellent extragalactic distance in- 
dicators. The quantitative analysis of their spectra - even at only moderate resolution 
- allows the determination of stellar parameters, stellar wind properties and chemical 
composition with remarkable precision. In addition, since the spectral analysis yields in- 
trinsic energy distributions over the whole spectrum from the UV to the IR, multi-colour 
photometry can be used to determine reddening, extinction laws and extinction. This is 
a great advantage over classical distance indicators, for which only limited photometric 
information is available, when observed outside the Local Group. Spectroscopy also al- 
lows to deal with the effects of crowding and multiplicity, as blue supergiants, due to 
their enormous brightness, are less affected by such problems than for instance Cepheids, 
which are fainter. 

Applying the FGLR method on objects brighter than My — —8 mag and using multi- 
object spectrographs at 8 to lOm-class telescopes, which allow for quantitative spec- 
troscopy down to my = 22 mag, we estimate that with 20 objects per galaxy we will be 
able to determine distances out to distance moduli of m — M ^ 30 mag with an accuracy 
of 0.1 mag. We emphasize that these distances will not be affected by uncertainties in 
extinction and metallicity, because we will be able to derive the corresponding quantities 
from the spectrum. 



7. Winds of hot massive stars 



All hot massive stars have winds, which are driven by radiation. As emphasized in 
the introduction these winds are fundamentally important for the spectral diagnostics 
of massive stars, for their evolution and for the galactic environment. Over the last 
decades very detailed and refined methods have been developed for the diagnostic of 
stellar winds and for modelling their hydrodynamics. A comprehensive review describing 
these methods and summarizing the basic properties of stellar winds was published a few 
years ago by Kudritzki & Puis (2000) (see also Kudritzki (2000) which was published 
in these Symposiiun Series). The mass-loss rates and terminal velocities of these winds 
are related to the physical parameters of massive hot stars through simple relationships. 
The stellar wind momentum Mvoo is related to radius R, luminosity L and metallicity Z 
through 



Mv^R^/^ oc Li-«(Z/Z0)-"-« (7.1) 
This the wind momentum luminosity relationship (WLR), which has been introduced 



first by Kudritzki, Lennon & Puis (1995) (Note that the left hand side is called the 
"modified stellar wind momentum" ). The physics leading to this relationship has been de- 
scribed in detail in Kudritzki (1998) and Kudritzki (2000) Proportionality constants for 
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for 0-stars and 



the WLR have been determined by |Puls, Kudrit zki, Herrero et al. (1996) 
by [Kudritzki et al. (1999), for B and A supergiants. (See also [Kudritzki fc Puis (20d0)| . 

The terminal velocities ■^oo of the winds of hot massive stars are related to the photo- 
spheric escape velocities vP^°* and metallicity through 



i;oo(x<:f (Z/Zo)-"-i5 (7.2) 

vP!!°' = (2GM(1 - r)i?-i)i/2 (7.3) 

The proportionality constant in Eg. 17.21 is about 2.6 (±0.5 ) for 0-stars and early B- 
supcrgiants but becomes smaller for lower temperatures and is about 1.0 for A-supergiants 
(see ,Kudritzki &: Puis (2000), but see also the new results discussed below). T is the usual 
distance to the Eddington limit. 

7.1. 0-star wind momenta based on new diagnostics with metal line blanketed model 

atmospheres 

The proportionality constants for the WLR were based on spectral diagnostics of Ha, 
which is very sensitive to the strengths of stellar winds and, therefore, regarded as a very 



good indicator of mass-los rates (see Kudritzki & Puis (2000)1. However, the original 
work was based on the use of NLTE model atmospheres, which neglected metal line 
blanketing. It was, therefore, crucially important to repeat the Hq studies of stellar 
winds with the improved model atmospheres described in the sections before. 



Repolust et al. (2004) and Markova et al. (2004) carried out a comprehensive study 
of 0-stars in the Milky Way. Fig. |S| indicates the general trend with regard to the di- 
agnostics of wind momenta. Since the analysis with metal line blanketed NLTE atmo- 
spheres yield lower temperatures and, therefore, lower luminosities, the WLR is simply 
shifted towards lower luminosities. Fig. E| shows the observed WLR for 0-stars. As in 
Kudritzki & Puis (2000) we include the results of similar diagnostics of Central Stars 
of Planetary Nebulae (CSPN) (for details see Kudritzki, Urbaneja, & Puis (2006a) I to 
demonstrate that very obviously the WLR extends to these low mass hot stars in a post- 



AGB evolutionary stage as well. Repolust et al. (2004) and Markova et al. (2004) discuss 
the difference between the WLR of dwarfs and supergiants and argue that the theory 
cannot reproduce this difference. They conclude that the Hq mass-loss rates of super- 
giants are affected by inhomogeneous stellar wind clumping and in reality are probably 
close to the ones of dwarfs. To check this hypothesis, we compare with the theoretical 
wind momenta for 0-stars and CSPN calculated by Kudritzki (2002) who used a new 



theoretical approach for the theory of line driven winds. The result is also shown in 
Fig. ^1 While we find a clear offset between the WLRs of dwarfs and supergiants, the 
theoretical wind momenta are a little too small in both cases, which would mean that 
clumping might be important for both dwarfs and supergiants. We will discuss the effects 
of clumping at the end of this review. 

With the new stellar parameters derived with metal line blanketed models there is also 
a change in the values of photospheric escape velocities. After a first check of all the new 
results the factor relating terminal velocity with escape velocity appears to be 3.1 rather 
than 2.6. 

7.2. Metallicity dependence 

Since the mechanism of driving stellar winds is the absorption of photospheric pho- 
ton momentum through many thousands of spectral lines, it is clear that the strengths 
of winds is expected to depend on metallicity. First predictions of the metallicity de- 



pendence of both mass-loss rates and terminal velocities were made by Abbott (1982) 
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Figure 16. Left: Observed wind momenta of O-stars and CSPN as a function of luminosity. 
Right: Calculated stell ar wind momenta for O-stars and CSPN using the theory of line driven 
winds as developed by [Kudri tzki (2002}] Symbols refer to model calculations. The dashed lines 
are the regression curves obtained from the observations displayed in the left diagram. The solid 
line represents the theoretical approach by [Vink et al. (2000)^ 




log(L/Ls,„) 



Figure 17. Left: Observed wind momenta of O-stars and CSPN as a function of luminos- 
ity. The shaded area represent observational results in the MW, LMC, and SMC and the 
solid lines are the observed WLR regression curves. The dashed curves represent theoretical 
predictions by 'Vink ct al. (2001)| for MW, LMC and SMC metallicities, respectively. (From 
IMokiem et al. (2006b)). Right: Calculated stell ar wind mom enta for O-stars and CSPN using 
tHc theory of line driven winds as developed by [Kudritzki (20 02) Symbols refer to model cal- 
culations. The solid curve is the observed regression for the S MC from the left d iagram. The 
dashed line represents the theoretical approach for the SMC by[Vink et al. (2001)[ 



and Kudritzki et aL (1987)| later confirmed and extended by Leitherer et al. (1992) 



Im- 



proved calculations were carried out more recently by Vink et al. (2001) and Kudritzki (2002) 
Observationally, [Puis, Kudritzki, Herrero et al. (1996) analyzed O-stars in the Clouds 
relative to the Milky way and found a first observational indication for the power law 
dependence of wind momenta on metallicity. Moreover, Kudritzki & Puis (2000) showed 
that O-stars in the metal poor SMC have lower lower terminal velocities than their MW 



counterparts. More recently, Massey, Bresolin, Kudritzki et al. (2004) 
Evans et al. (2004)[ [Hillier et al. (2003)[ and [Crowther et al. (2002)[ have clearly con- 
firmed these results. l''ig.ll7lshows the most recent results obtained by [Mokiem et al. (2006b)[ 
which indicate excellent agreement with the predictions of the theory. 



Massey, Puis, Pauldrach et al. (2005)] 
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Figure 18. Left: Observed wind momenta of Bla supergiants as a function of luminosity for 
objects hotter than 23kK (BO-BO. 5) (upper panel) and cooler than 23kK (BO. 7-3) (lower panel) 
and for objects in the Galaxy (circles), LMC (squares) and SMC (triangles). Overplotted are wind 
momenta predicted b y the computations by Vink et al. (20 00, 2001) for the metallicities of the 
three galaxies, (from jCrowther, Lennon fc Walborn(2005)| l. Right: Mass- loss rates as predicted 
by Vink et al. (2001) for SMC B-supergiants compared to observations. Circles represent early 
type supergiants, stars corresponds to mid types. 



7.3. B supergiants 

Because of changes in ionization leading to different sets of spectral lines absorbing 
photospheric photon momentum and driving the stellar winds the WLR is expected to 
be spectral type dependent. Indeed, Kudritzki et al. (1999) in their study of winds of B 
and A supergiants found a strong variation of the WLR with spectral type. While early 
B supergiants seemed to have wind momenta only somewhat weaker than their 0-stars 
counterparts, mid B supergiants showed much weaker winds, whereas A supergiants had 
stellar wind momenta comparable to the early B spectral types, but a steeper slope of 

the WLR. 

With the new line blanketed models available Trundle et al. (2004)[ Trundle &: Lennon (2005)1 
Urbaneja (2004)[ Evans et al. (2004) and Crowther, Lennon fc Walborn(2005) have re- 



investigated the wind properties of B-supergiants. In general, the results by Kudritzki et al. (1999) 



are confirmed. There is still on average a difference in the strengths of wind momenta 
between early and mid B-spectral types. However, while mass-loss rates for early spectral 
types are very similar to those obtained with the unblanketed models by Kudritzki et al. (1999)[ 
the new values for mid B-supergiants are about a factor of 3 higher on average, but even 
with this increase they remain significantly lower than those for the earl B types (see 
Fig. I18|) . There is a clear metallicity dependence of wind momenta. 

As shown in Fig.^lthe theory of line driven winds fails to reproduce these observations. 
Vink et al. (2000, 2001) predict strong increase of mass-loss rates below 24kK due to 
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changes of ionization, which is definitely not observed. For the early type supergiants 
show mass-loss rates much stronger than predicted by the theory. While this can be 
explained by stellar wind clumping, we believe that that in general more theoretical 
work for B-supergiant winds is needed. 

The ratio between terminal and escape velocities changes as a function of effective 
temperature showing a similar trend as found by Prinja fc Massa (1998)| but quantita- 
tively different. For temperatures above 24kK the ratio is 3.4, between 20kK and 24kK 
it is 2.5 and below 20kK 1.9 is found. In all temperature ranges there is a large scatter 
around these average values, though. 

7.4. Winds at very low metallicity 

As indicated in the introduction and as discussed in several contributions throughout 
this symposium there is growing evidence that the evolution of galaxies in the early 
universe is heavily influenced by the formation of first generations of very massive stars. 
Thus, it is important to understand the nature of radiation driven winds at very low 
metallicity. As has been shown by Kudritzki (2002)[ radiative line forces show a different 
dependence on optical depth and electron density at very low metallicity, which requires 
significant modfications of the theoretical description. With those implemented it can be 
shown that the power law dependence on metallicity of mass-loss rates, wind momenta 
and velocities breaks down and a much stronger dependence on metallicity is found. For 



details we refer the reader to the paper by Kudritzki (2002) which also provides ionizing 
fluxes and predicted UV spectra as a function of metallicity. 

7.5. The effects of rotation and instabilities 

There are many mechanisms, which might lead to an enhancement of mass-loss during 
stellar evolution. A classical mechanism already taken into account by |Pauldrach, Puis, fc Kudritzki] 
and Friend & Abbott (1986) is through the centrifugal forces provided by stellar rotation. 
Kudritzki fc Puis (2000)1 [Petrenz fc Puls(2000)| 



and 



Owocki (2005) 



give 



an overview 



about the possible effects. In particular at low metallicities, when massive 0-stars are hot- 
ter and have much smaller radii and the effects of radiation driven winds become smaller, 
stellar rotation might become a crucial mechanism not only for stellar mass-loss, but also 
for stellar evolution (see|Marigo eFal. (2003) , Maeder et al. (2005)||Meynet et al. (2005)| 



Hirschi et al. (2005)|pappim et al. (2006)|jMeynet et al. (2006)f . 



Continuum driven winds and their instabilities as very likely encountered in LBV out- 
bursts are another important physically mechanism at least for objects very close to the 
Eddington limit (seepiaviv (2001)[|Shaviv (2005)] [Owocki et al. (2004)[|Owocki (2005)| 



Smith fc Owocki (2006) and the contribution by Nathan Smith at these proceedings). 
Stellar pulsations have also been discussed frequently to enhance stellar mass-loss, but 



as shown by Baraffe et al. (2001) they become less important at low metallicities. 



7.6. The problem of the "weak wind stars" 
UV and optical studies of 0-stars in the SMC with luminosities L < 10^-^ Lq by Bouret et al. (2003) 



and Martins et al. (2004) indicated much smaller stellar wind momenta than expected 
from the theory of radiation driven winds and from a simple extrapolation of the WLR 
from higher luminosities. More recently, Martins et al. (2005b) have investigated similar 
"weak wind stars" in the Milky Way and it was found that 0-dwarfs with L ^ 
indeed seem to have much weaker winds than predicted by the theory with a discrepancy 
of the order of a factor of hundred. It is always simple to speculate about too strong 
stellar winds. One can invent additional wind driving mechanisms or blame the neglect 
of clumping in the diagnostics as the reason. However, it is very difficult to explain winds 
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Wavelength [A] Teff [kK] 

Figure 19. Left: Radiative transfer fits of the far UV PV lines of MW 0-stars. Right: Ionization 
fraction of PV as deter mined from fits of Ha and PV as function of effective temperature. (From 
[FuUcrton ct ah (2006) 



that are too weak, because the radiative force at a given luminosity is always there and 
cannot be simply switched off. Thus, it is natural to be suspicious about the accuracy of 
the spectroscopic diagnostics of these objects. 

The determination of mass-loss rates comes mostly from the analysis of UV metal 
lines, since Hq provides only upper limits in most of the cases. While the diagnostics of 
those lines have been done in the most careful way with the state-of-the-art model atmo- 
spheres described before, the mass-loss rates determined depend crucially on ionization 
calculations which could severely be affected by the soft X-ray emission of shocks in the 
stellar wind flow. The authors are aware of this problem and have included effects of 
shock emission, but it is open at this point whether this treatment of shocks is sufficient. 
Other possible options - from our point of view less likely - to explain the weak wind 



stars are discussed by Martins et al. (2005b) and Mokiem et al. (2006b) 



7.7. Stellar wind clumping 
While Hq is, in principle, a perfect tool to measure mass-loss rates (see Kudritzki fc Puis (2000)| 



for discussion and references), the results might be affected by stellar wind clumping. It 
has been known since long that line driven winds are intrinsically unstable (Owocki et 
al. 1988, 2004). This might lead to inhomogeneous, clumped winds such as described 



by Owocki & Runacres (2002) with regions of enhanced density pd and regions, where 
the density is much lower. In a very simple description, introducing clumping factors fd 
similar as in PN diagnostics, the relationship between the average density of the stellar 
wind flow pav and the density in the clumps is then given by pd = Pavfci- The same 
relationship holds for the occupation numbers Ui of ions. 

Line opacities n depend on density through k cc Hi (x p^ and for very small, optically 
thin clumps the avarage optical line depth in the wind is given by Tav oc nf"" oc nf f^^ oc 
Pavf^~^- For a dominating ionization stage we have x — 1 and the clumping along the 
line of sight cancels and does not affect the diagnostics. However, bound hydrogen is a 
minor ionization stage in hot stars depending on recombination from ionized hydrogen 
with ni{H) oc nsnp oc p^. Thus, if fd is signiflcantly larger than one, the Hq, mass-loss 

1/2 

rate diagnostic is systematically affected and we have M{Ha) — M{true)fJ , following 
from the fact that M{true) oc pav 

The spectral diagnostics of clumping is difficult. In principle, it requires the compar- 
ison of lines with different exponents x in the density dependence of their opacities. In 
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Figure 20. Diagnostic of clumping factors in the winds of CSPN. The first two rows show model 
fits of the (left) and Hell 4686 (right) profile of IC418 with f^i = 1 (top row) and fd = 50 indi- 
cating enormous clumping. The bottom row shows the same for He2-108 and fd = 1 demonstrat- 
ing that for this object the wind is very likely homogene ous. Note that the Hg and Nil nebula r 
emission lines have been (imperfectly) subtracted. (From|Kudritzki, Urbaneja, fc Puis (2006a)|l. 



WR-type stars with very dense winds and very strong wind emission lines incoherent 
electron scattering produces wide emission wings, the strength of which goes with x ^ I. 
Clumping factors of the order of ten to twenty were found ( |Hillier (1991)1 see also contri- 
bution by Paul Crowther in these proceedings). This technique does not work for 0-type 
stars, as their winds have much lower density. Also the UV P-Cygni lines of dominating 
ions provide usually little help, as these lines are mostly saturated and the ionization 
equlibria are uncertain. However, in most recent work on massive 0-stars using FUSE 
and Copernicus spectra the PV resonance line at 1118 and 1128 A has been used as 
an indicator of clumping. The advantage of PV is the low cosmic abundance so that 
the line is completely unsaturated even when in a dominating ionization stage. Substan- 



tial clumping was found ( Hillier et al. (2003)| Bouret et al, 2003, 2005) with clumping 
factors of the order of ten. 

Very recently, 'FuUerton et al. (2006) have carried out a comprehensive study of Milky 
Way 0-stars with well observed FUSE PV line profiles producing detailed PV radiative 
transfer line fits to determine the product q{P*^)M, where g(P'*^) is the ionization 
fraction of PV in the ground-state (see Fig ll9|) . Comparing with mass-loss rates derived 
from Hq, (or radio free- free emission) they produced the plot of g(P'*+) as function of 
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Te// also shown in Fig^l Assuming that PV is a dominating ionization stage at all 
temperatures (and, therefore, g(P'*+) = 1), they conclude that the Hq mass-loss rates 
are too high by a very large factor up to many orders of magnitudes and implying 
enormous filling factors. 

This is very important work based on the best available diagnostic techniques and poses 
a very serious problem. The crucial assumption is, of course, g(-P^^) = 1. Test calcula- 
tions done by us with the model atmosphere code FASTWIND show that models (without 
shock emission) predict g(P*+) — 1 only for Tg^^ ^ ibkK, which would imply a much 
lower effect then at hotter temperatures. Clearly, for future work a very detailed investi- 
gation of the PV ionization (including shock emission) is needed to address this funda- 



mental problem of clumping in 0-star winds. The results found by FuUerton et al. (2006) 
are certainly alarming. 

There is an alternative method for the diagnostics of clumping at least for cool 0-stars 
with Te// ^ 37, OOOK, where Hell is a dominant ionization stage. That means for objects 
with strong winds and Hell 4686 in emission and formed in the wind this line should 
have a density dependence close to x = 1. Its relative strenght to Hq, should allow to 
constrain id- 



Kudritzki, Urbaneja, & Puis (2006a) have applied this technique to study clumping in 
the winds of CSPN with very interesting results (see Fig. I20() yielding clumping factors 
varying in a range from 50 to 1. 
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